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Abstract. We present component-resolved and total pair distribution functions for 
a 2DEG with two symmetric valleys. Our results are based on quantum Monte Carlo 
simulations performed at several densities and spin polarizations. 

O . PACS numbers: 71.10.Ca, 71.45. Gm, 02.70.Ss 

o 



1. Introduction 

A two-valley (2V) two-dimensional (2D) electron gas (EG) is the simplest model to 
describe electrons confined in solid state devices such as Si-MOSFETs [T] and certain 
quantum wells [2]. In this 2D model, electrons interact via a 1/r potential in a uniform 



neutralizing background and possess an additional discrete degree of freedom (the valley 
index). One may identify electrons with given spin projection and valley index as 
belonging to a different species or component. We focus on the case of two symmetric 
valleys, where the number of up (down) spin electrons is the same for both valleys. 



In this case at zero temperature the 2DEG is completely characterized in terms of the 
coupling parameter r s = 1/ 'y/n na^ and the spin polarization ( = (n-f — n>i)/n (where 
n is the total electron density, ob the Bohr radius, n^Q) the density of up (down) spin 
electrons) . 

The interest in the properties of the 2DEG has been strongly revived in the last years 
due to the experimental discovery of a previously unexpected metal-insulator transition 
[31 H] in which the valley degree of freedom appears to play an important role [5j [6]. The 
transition takes place at low density, where an accurate treatment of electron correlation 
is crucial. In this respect, quantum Monte Carlo (QMC) simulations have provided over 
the years the method of choice for microscopic studies [3 [9], [101 [11], [121 OS] °f the 
2DEG. 

We have recently provided an analytic expression of the correlation energy of the 
2V2DEG [H]. Here, we focus on the pair distribution functions (PDFs), which are 
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strictly related to the description of exchange and correlation properties of the system 
(see e.g. Ref. [15]). As pointed out in Ref. [13], PDFs may serve a variety of purposes, 
among which the estimate of finite thickness effects on the 2DEG spin susceptibility 
[16], applications in DFT calculations, or a test of the accuracy of hypernetted-chain 
calculations [T7| . 

2. Pair distribution function: definitions 

The PDF is related to the probability of finding two electrons at positions r and r' 
respectively. The component-resolved PDF g a/ g(r',r) of a multicomponent 2DEG is 
defined as [33] 



with if)^, ip a denoting creation and annihilation field operators, (. . .) the expectation 
value on the ground state and n a (r) = (^(r)^ a (r)) the electron density of the 
component a. The normalization is such that g a p = 1, in case there is neither exchange 
nor correlation. In a homogeneous and isotropic system, g a p depends only on the relative 
distance r = |r — r'| and there is symmetry for index permutations (g a p = gp a )- If 
c a = n Q /n is the concentration of the component a, the total (component- summed) 
PDF g(r) reads 



In the 2V2DEG, a = ov is a composite index which denotes the spin and valley 
(respectively a and v) degrees of freedom and spans the four cases a =| 1, \ 1, | 2, \ 2. 
In general, for 2V there are ten different g a p, but in the case of two symmetric valleys 
the number of different g a p is two for ( = 0, 1 and five for < ( < 1 (see also Ref. [T7]). 
In the following we shall denote the different g a p with one of the possible different labels 
(e.g. for < C < l gnu = 9m? = 9ni2 = 9mi)- 

3. Simulation details 

Most of the simulation details are the same as in Ref. [14] . We performed fixed-phase 
DMC simulations (for a review of QMC techniques see e.g. Ref [18]) for r s = 1, 2, 5, 10, 20 
and ( = 0,3/13,5/13,7/13,10/13,1. To reduce the finite size effects, we used twist- 
averaged boundary conditions (TABC) [TJJ], which also allow to change ( by flipping 
any number of spins at fixed number of electrons N. To sample the PDFs we performed 
simulations for a system of N = 52 electrons. Time steps were chosen at the different 
r s to give an acceptance rate corresponding to ~ 99%. We did simulations with 320 
walkers. The twist grid is the same as in Ref. |14j . As in Ref. [14J, we used a Slater- 
Jastrow trial wave function but, here, we considered only plane-wave nodes, since 
the more accurate backflow (BF) nodes yield only slight modifications to the PDFs 





(1) 



n (r')n a (r) 



9( r ) =J2 c <* c P9ap(r). 
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[8], 02]. Besides, BF effects on the energy were shown to be bigger in the two- component 
case than in the four- component system |14j . 

DMC provides the mixed estimate of an operator O, i.e. (0) m i X = 
(^o\0\^t)/ (^o^t) (with \l/ denoting the ground state of the system). If O commutes 
with the Hamiltonian H, (0) m i X coincides with the expectation value on the true ground 
state \EV If O does not commute with H (as in the case here considered), it is better to 
compute the extrapolated estimate (0) extr = 2(0) mix — (0)vmc + 0(5 2 ) (where (O)vmc 
is the variational MC expectation value on *&t) which is accurate to second order in the 
difference 6 between ty and *&t- 

4. Two valley pair distribution functions 




In Fig. [I]|2] we show representative examples of 2V component-resolved PDFs for 
C = 0, 1 and finite ( respectively. All lenghts are given in units of a = r s aB- The 
component-resolved PDFs shown in Fig. [1] illustrate the tendency to a local order which 
favors electrons belonging to different species (as e.g. f 1 and j 1 ) to get closer than 
electrons belonging to the same species. For intermediate spin polarizations (see Fig. [2]) 
the component-resolved PDFs exhibit a richer structure than the ( = 0, 1 cases, with 
qualitative features clearly related to the interplay of exchange and correlations (for 
example, in the density range considered the diagonal PDF of the minority component in 
a strongly polarized system is found to be determined by exchange alone, to a very good 
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Figure 2. Examples of 2V component resolved pair distribution functions at finite (. 
See labels. 



approximation). The significant spin-polarization dependence seen in the component- 
resolved PDFs almost disappears in the total PDFs, particularly at large r s . This can 
be appreciated in Fig. [3] (left panel), which shows the total PDFs for a high density 
(r s = 1) and a low density (r s = 20) case with zero and full spin polarization. The 
dependence of the total PDF on r s is depicted in Fig. [3] for ( = 0.538. 

We see that even for r s as small as 1 the effect of exchange on the total PDF, 
which is expected to decrease with the number N c of (equally populated) components, 
is already very small for N c = 2. 




0123456 0123456 



r/a r/a 

Figure 3. 2V total pair distribution functions. Left panel: r s = 1,20 and C = 0, 1. 
Right panel: r s — 1, 2, 5, 10, 20 and £ = 0.538 (increasing peaks for increasing r s ). 
All lengths are in units of a = r s aB- 
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Full tabulations of the calculated PDFs are available upon request from the first 
author. 
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